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^ — , ■ In this paper we study isoperimetric problems of the calculus of variations with left 

and right Riemann-Liouville fractional derivatives. Both situations when the lower 
bound of the variational integrals coincide and do not coincide with the lower bound 
\0 ■ of the fractional derivatives are considered. 
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1 Introduction 



Isoperimetric problems consist in maximizing or minimizing a cost functional subject to in- 
tegral constraints [5j. They have found a broad class of important applications throughout 
the centuries. Areas of application include astronomy, physics, geometry, algebra, and anal- 
ysis [6l [T7j . Concrete isoperimetric problems in engineering have been also investigated by a 
number of authors |18| . 

The study of isoperimetric problems is nowadays done, in an elegant and rigorously way, 
by means of the theory of the calculus of variations. This is possible through a powerful 
tool known as the Euler-Lagrange equation |33| . Recently the theory of the calculus of 
variations has been considered in the fractional context 0[51[ini[IIl[I21[i[I2E3[23[271[33]. 
The fractional calculus allows to generalize the ordinary differentiation and integration to an 
arbitrary (non-integer) order, and provides a powerful tool for modeling and solving various 
problems in science and engineering [28 1 129 1 131] ■ The problems considered are more general, 
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and hold for a bigger class of admissible functions which are not necessarily differentiable in 
the classical sense [30]. Several results were proved for the new calculus of variations. They 
include: Euler-Lagrange equations for fractional variational problems with Riemann-Liouville 
PP, Riesz [3], Caputo, and (a, (3) derivatives [19]; transversality conditions [2]; and Noether's 
symmetry theorem [21\ I22|. For a state of the art of the fractional variational theory see the 
recent papers [H [151 EH EH 121] and references therein. In this paper we develop further 
the theory of the fractional variational calculus by studying isoperimetric problems. 

The paper is organized as follows. In Section[2]we shortly review the necessary background 
on fractional calculus. Our results are given in Section[3j In Section f3. II we introduce the basic 
fractional isoperimetric problem and prove correspondent necessary optimality conditions, 
both for normal and abnormal extremizers (Theorems [6] and respectively). In Section [3.21 
we generalize our results for functionals where the lower bound of the integral is greater than 
the lower bound of the Riemann-Liouville derivatives. Finally, in Section 13.31 we present a 
necessary condition of optimality for the case where the order of the derivative is taken as a 
free variable. 

2 Preliminaries of fractional calculus 

A fractional derivative is a generalization of the ordinary differentiation, which allows real 
number powers of the differential operator. There exist numerous applications of fractional 
derivatives to several fields, like geometry, physics, engineering, etc. In the literature we 
may find a great number of definitions for fractional derivatives (see, e.g., [281 \29\ I31j). In 
this paper we deal with the left and right Riemann-Liouville fractional derivatives, which are 
defined in the following way. 

Definition 1. Let f : [a,b] — >• K. be a continuous function. The left and right Riemann- 
Liouville fractional derivatives of order a > are defined respectively by 

aV«f(x) = - — / (x - t) n - a - l f(t)dt , x £ (a, b] , 

1 [n — a) ax n J a 

and 

x V%f(x)= ) > — / (t-x) n -«- l f(t)dt, xe[a,b), 
1 (n — a) ax n J x 

where T is the Euler gamma function, a is the order of the derivative, and n = [a] + 1 with 
[a] being the integer part of a. 

If a > 1 is an integer, these fractional derivatives are understood in the sense of usual 
differentiation, that is, 

cP2f[x) = (£) f^) and x Vtf{x) = L±Y f( x ). 

^From the physical point of view, if f(x) describes a certain process through time x, then 
the left derivative is related to the past of this process, while the right derivative belongs to 
the future. 

These operations are linear, in the sense that 

a V^f(x) + ug{x)) = n a V«f{x) + v a V«g{x) 
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and 

x VZ(vf{x) + vg(x)) = ti x V%f(x) + v x Vtg{x). 
We now present the integration by parts formula for fractional derivatives. 

Lemma 2. ( |314 p. 46]) If f and g and the fractional derivatives oX^xd an d x^bf are contin- 
uous at every point x £ [a, b], then for < a < 1 we have 

b ,-b 
f(x) a V«g(x)dx = / g(x) x VZf(x)dx. (1) 



Moreover, formula fi]J is still valid for a = 1 provided f or g are zero at x = a and x = b. 

3 Main results 

^From now on we fix a, f3 € (0, 1). We consider functionals J of the form 

J(y)= / L(x,y, a V%y, x VPy)dx (2) 

J a 

defined on the set of admissible functions y that have continuous left fractional derivatives of 
order a and continuous right fractional derivatives of order j3 in [a, b], and where (x, y, u, v) — > 
L(x,y,u,v) is a function with continuous first and second partial derivatives with respect to 
all its arguments such that y, aDxVi x^^y) has continuous right fractional derivative of 

order a for all x € [a, b] and ^ (x, y, aVxVi xT^bV) nas continuous left fractional derivative of 
order f3 in [a, b\. 

Remark 1. The left Riemann-Liouville fractional derivative is infinite at x = a if y(a) ^ 0. 
If y{°) 7^ 0; then the right Riemann-Liouville fractional derivative is also not finite at x = b 
[30]. For this reason, by considering that the admissible functions y have continuous left 
fractional derivatives, then necessarily y(a) = 0; by considering that the admissible functions 
y have continuous right fractional derivatives, then necessarily y(b) = 0. This fact seems 
to have been neglected in some previous work on the calculus of variations with Riemann- 
Liouville fractional derivatives. Alternatively, we can consider the general case of boundary 
conditions, say y{a) = y a and y(b) = y\,, and study functionals of type 

rb 

J(y) = / L (x, y{x), oPxivi.x) - y a ), x V?(y(x) - y b ))dx. 



b 
J a 

This needs, however, a modified fractional calculus \26§ . 

Definition 3. The functional J is said to have a local minimum (resp. local maximum,) at 
y if there exists a S > such that J{y) < J{y\) (resp. J{y) > J{y\)) for all y\ such that 
\\y - 2/1 II < S. 

In [lj the following problem is addressed: among all curves y(x) satisfying the boundary 
conditions, find the ones that maximize or minimize a given functional J . An answer to this 
question is given in the next theorem. 

Theorem 4 (|lj). Let J be afunctional as in (dJJ and y an extremum of J . Then, y satisfies 
the following Euler- Lagrange equation: 

+ . sf _ + . BS _ = a (3 ) 



3 



3.1 The fractional isoperimetric problem 

We introduce the fractional isoperimetric problem as follows: find the functions y that satisfy 
boundary conditions 

y{a) = Va, y{b) = Vb (4) 

(Va = if left Riemann-Liouville fractional derivatives are present in ([2]); y^ = if right 
Riemann-Liouville fractional derivatives are present in ([2]) — cf . Remark [1]) , the integral 
constraint 



Ay) = \ 9{x,y, aPxV, x T%y)dx = I 

J a 



(5) 



and give a minimum or a maximum to ([2]). We assume that I is a specified real con- 
stant, functions y have continuous left and right fractional derivatives (if present in ©), 
and (x, y, u, v) — > g(x, y, u, v) is a function with continuous first and second partial deriva- 
tives with respect to all its arguments such that ^{x,y, oVxVi xD^y) has continuous right 
fractional derivative of order a for all x G [a, b] and aVxVi x^y) nas continuous left 

fractional derivative of order j3 in [a, b\. Theorem |4] motivates the following definition. 

Definition 5. An admissible function y is an extremal fori in if it satisfies the equation 



dy +x b du +a x dv 



for all x £ [a, b] . 



The following theorem gives a necessary condition for y to be a solution of the fractional 
isoperimetric problem defined by ([2])-(j3|)-([5]) under the assumption that y is not an extremal 
for 1. 

Theorem 6. Suppose that J given by (0j has a local minimum or a local maximum at y 
subject to the boundary conditions |J) and the isoperimetric constraint {5|j. Further, suppose 
that y is not an extremal for the functional X. Then there exists a constant A such that y 
satisfies the fractional differential equation 

dy + xVb du + aV * dv ~ ° (6) 

with F = L — Xg. 

Proof. Consider neighboring functions of the form 

y = y + ei7/i + e 2 r/ 2 , (7) 

where for each i G {1,2} e, is a sufficiently small parameter, r/j have continuous left and right 
fractional derivatives, and r]i(a) = r]i(b) = 0. 

First we will show that ([7]) has a subset of admissible functions for the fractional isoperi- 
metric problem. Consider the quantity 



,-b 

I g(x, y + exrjx + e 2 r]2,aVxV + e ia^xVi + ^aPx^x^V + ^lx^bVi + ^x^^dx. 

J a 
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Then we can regard T{y) as a function of ei and e 2 . Define J(ei, e 2 ) = I(y) — I- Thus, 



7(0,0) = 0. 



(8) 



On the other hand, we have 



de 2 



(0,0) 



b \dg dg 



dg_ 

dv 

8y +x b du +a x dv 



>l2 + Z» P "'fe +j^aP$m 



•q 2 dx, 



dx 



(9) 



where ([9]) follows from ([I]). Since y is not an extremal for X, by the fundamental lemma of 
the calculus of variations (see, e.g., [33l p. 32]), there exists a function 772 such that 



dl_ 

de 2 



^0. 



(10) 



(0,0) 



Using ([8]) and (|10p . the implicit function theorem asserts that there exists a function e 2 (-), 
defined in a neighborhood of zero, such that j(ei, e2(ei)) = 0. We are now in a position to 
derive the necessary condition ([6]). Consider the real function J(ei, e 2 ) = J{y)- By hypothesis, 
J has minimum (or maximum) at (0, 0) subject to the constraint 1(0, 0) = 0, and we have 
proved that V 7(0, 0) 7^ 0. We can appeal to the Lagrange multiplier rule (see, e.g., [33, 
p. 77]) to assert the existence of a number A such that V( J(0, 0) — AJ(0, 0)) = 0. Repeating 
the calculations as before, 



and 



Therefore, one has 



8J_ 

de~l 



(0,0) 



(0,0) 



dL dL a dL 

8y~ +x b d^ + a x dv~ 



dy +x b 8u +a x dv 



rj\(x)dx 



rj\{x)dx . 



dL 



.dL 



dL 



A 



dy b du a ~ x dv \dy 
Since (jlip holds for any function 771, we obtain 



+ X V b du +aV XQv 



rfr(x)dx = 0. (11) 



dL 



dL 



; dL 



— + X V%— + a Vi\— - A + + a V^ 



dy 



du 



dv 



idg 



dy 



du 



dv 



□ 



Remark 2. Theorem [6] holds true in the case when a or (3 are equal to 1. Indeed, in the 
proof we imposed the condition 772(a) = 772(b) = 0, and formula fTP is valid. 
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Example 1. Let a be a given number in the interval (0, 1) 
fractional isoperimetric problem: 



We consider the following 



[\x 4 + ( V%y) 2 )dx 
Jo 



mm 



x 2 T>"y dx = - 
o 5 



2/(0) = 0, y(l) 



2a + 3a 2 + a 3 ' 
The augmented Lagrangian is 

F(x, y, V«y, x v{y) = x 4 + ( V«y) 2 - A x 2 V«y 

and it is a simple exercise to see that 



y{x) 



1 



-dt 



1 



2x 



a+2 



r(a) J (x - t) 1 -* r(a) 2a + 3a 2 + a 3 



(12) 



(13) 



(i) is not an extremal for the isoperimetric functional, (ii) satisfy oD^y = x 2 , (Hi) (0) holds for 
A = 2, i.e., x T>i (2 oTJ x y — 2x 2 ) = 0. We remark that for a = 1 (\13h gives y(x) = x 3 /3, which 
coincides with the solution of the associated classical variational problem (Fig.\7ty. Indeed, for 






= 0.2 




= 0.5 


/ 




/ 




/ 

,a 


= 0.7 


a 


= 0.9 


//Ct 


= 1 



X 



0.25 0.5 0.75 1.0 

Figure 1: The fractional solution converges to the classical one as a — > 1. 



a — y 1 our fractional problem tends to the classical isoperimetric problem of minimizing 
the functional f (x 4 + {y') 2 )dx subject to the isoperimetric constraint Jq x 2 y' dx = | and the 
boundary conditions y(0) = and y(l) = 1/3. Then, F = x 4 + (y') 2 — Xx 2 y' and the classical 
Euler- Lagrange equation is 



dy dx (ch/ ^ ^ 



0. 



(14) 



The solution of J 14\ ) subject to y(0) = 0, y(l) = 1/3, and x 2 y' dx = | is A = 2 and 
y = x 3 /3. 
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Introducing a multiplier Ao associated with the cost functional ([2]), we can easily include 
in Theorem [6] the situation when the solution of the fractional isoperimetric problem defined 
by ([21)-©-© is an extremal for the fractional isoperimetric functional. This is done in 
Theorem [7J 

Theorem 7. If y is a local minimizer or a local maximizer of (QJ) subject to the boundary 
conditions ^ and the isoperimetric constraint |5]), then there exist two constants Xq and X, 
not both zero, such that 

OK n dK R dK , . 

■"far-'-* 1 **- < 15 > 



Oy 



with K = XqL — Xg. 

Proof. Using the same notation as in the proof of TheoremEl we have that (0, 0) is an extremal 
of J subject to the constraint 1 = 0. Then, by the abnormal Lagrange multiplier rule (see, e.g., 
[33| p. 82]) there exist two reals Ao and A, not both zero, such that V(AoJ(0, 0) — A/(0, 0)) = 0. 
Therefore, 



Ar 



dJ_ 

de 1 



X 



(0,0) 



di_ 

de 1 



0. 



(0,0) 



Applying the same reasoning as in the proof of Theorem El we end up with ([15 



□ 



3.2 An extension 

In [9] a fractional functional 



C(y) 



L(x,y, a V%y)dx 



(16) 



is considered with [A, B] C [a, b], i.e., with the lower bound of the integral not coinciding with 
the lower bound of the fractional derivative. The main result of [9] is a new Euler-Lagrange 
equation for the functional (j!6p . We now extend the techniques of [9j to prove an Euler- 
Lagrange equation for functionals containing both left and right Riemann-Liouville fractional 
derivatives, i.e., for fractional functionals of the form 

J(y)= L(x,y, a V«y, x V?y)dx, (17) 
J A 

where the integrand L satisfies the same conditions as before. Let y be a local extremum of 
J such that y(a) = y a and y{b) = yt, and let y = y + erj with rj(a) = n(b) = 0. Consider the 
function J(e) = J{y + en). Since J(e) has a local extremum at e = 0, then 







dL dL _ dL a 



jV^ndx 



+ 



— dx+\ I — 
dy [Ja du 

b dL a , f b dL n : , 
— • x T>< h r}dx - — ■ x V^r]dx 
A 9v J B dv 



dx 

A dL 

du 



" B dL , 
— — • rjdx + 



J A dy 

+ I VAD p x -^dx 



B 



dL 



V-xV B Q^dx 



V ■ BUt-^dx 



aV^ndx 



dL 



11 ' xVA d^ dx 



B 



dv 
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Continuing in a similar way, 



f B dL 


Ja 


dy 






+ 


\f 




J A 




4 r 



rA qj^ rB gj^ i-A qj^ 



V ' a^x "ET~ + / V ' aT^x tt - dx 



dL 



dv 

pa dL dL 



B 



dv 



+ 



^dL p dL 
A ~dv ~ B ov 



rjdx + 
rjdx . 



B 



dL 
dy 



^edL , 
V ■ B^—dx 
B dv 



Let r/i : [a, A] — > R be any function satisfying 771(a) = 0, and r/ be given by 

r](x) = 



r]i(x) if x G [a, A] , 
elsewhere. 



Therefore, 







V 5L ^ 



rj\dx. 



By the arbitrariness of ryi and the fundamental lemma of calculus of variations, 

du x A du 



dl dl 
xV%^- - X V%— = for all x G [a, A]. 



Analogously, we have 



and 



dT dT dT 

-q^+ *V%^+ = for all x G [A, 5] , 



A^f ^ - bV^ = for all x G [S, 6]. 



<9t> 



dv 



We have just proved the following. 

Theorem 8. Let y be a local extremizer of fli?| ). Then, y satisfies the following equations: 

' dL dL dL 

■^ + * V b^+aVP— = forallxe[A,B], 

dL „ n dL 



xVB du 



dv 



dv 



for all x G [a, A] , 
for all x G [B,b]. 



Remark 3. Theorem^ simplifies to the result proved in f^j in the case the Lagrangian L in 
(17) does not depend on the right Riemann-Liouville fractional derivative xT^hV- 

We will study now the fractional isoperimetric problem for functionals of type ()1T|) subject 
to an integral constraint 



rB 

x {y) = j 9(x, y, a L>xV, xL>^y)dx = I . 
J A 



(18) 
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Definition 9. We say that y is an extremal for functional I given in \18)) if 

for all x € [A,B\. 



d 9 , ~ a d 9 , ^p d 9 
dy :+xVB d^ + AVx di 



Theorem 10. Let y give a local minimum or a local maximum to the fractional functional 
[11) subject to the constraint U8\) . If y is not an extremal for I , then there exists a constant 
A such that 

' BF BF BF 

-^+^B^+^ x — = forallxe[A,B] 

BF BF 

for all x € [a, A] (19) 



T) a 



Bv 



x Bv 



for all x € [B,b] 



with F = L — Xg. 



Proof. Consider a variation (€1,62) ^ y = y + + 62772 where 771(a) = 771(6) = 772(a) 
772(6) = 0. Let 



/(ei,e 2 ) 



g(x,y, aVxVi xV h i))dx - I. 



Then, 1(0, 0) = and 



Be 2 



(0,0) 



B 



B 9 . dg a Bg « 

^2 + ^ 772 + ^772 



Bg 
' B ~Bu~ 



x u B~ x u A 



Bg_ 

Bu 



+ 



B 



x Bv x Bv 



mdx + 
n 2 dx. 



dx 

B 'dg 
By 



Ov 



r]2 dx 



Since y is not an extremal for I, there exists a function 772 such that 7^ 0. By the 



Oc'2 



(0,0) 



implicit function theorem, there exists a subset of curves {y + 61771 + 62772 | (£1,^2) £ K 2 } 
admissible for the fractional isoperimetric problem. Let j(ei,62) = <J(y)- Then, there exists 
a real A such that V(J(0,0) - A/(0,0)) = 0. Because 



BJ_ 



BI_ 

dTi 



(0,0) 



v adL BL 

xVb b^~ xVa b^ 



+ 



B 



(0,0) 



Bg 

xUB -Bu~~ xL/A Bu 



+ 



B 



8J_ 

06! 



(0,0) 



rjidx + 
r/idx, 

r]idx + 

7/1 dx, 

BI_ 

de 1 



B 



BL 



BL 



BL 

x Bv 



r\\dx 



B 



99 



Bg 



b^ + xVb b^ + aVx b^ 



r]idx 



A 



0, 



(0,0) 



and 771 is an arbitrary function, it follows (|19|) . 



□ 
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Similarly as before, we can include in Theorem 1101 the situation when the solution y is an 
extremal for X (abnormal extremizer). For that we introduce a new multiplier Ao that will be 
zero when the solution y is an extremal for I and one otherwise. 

Theorem 11. If y is a local minimizer or a local maximizer of |i7| ) subject to the isoperimetric 
constraint \18\) . then there exist two constants Xq and A, not both zero, such that 

( clK BK r)K 

% + ^¥ + ^¥ = ° forallxe[A,B} 



dK 
du 



dK 



for all x E [a, A] 
for all x € [B, b] 



with K = XqL — \g. 



3.3 Dependence on a parameter 

Consider the following fractional problem of the calculus of variations: to extremize the 
functional 

" ' x a 



{0 D«y) 2 -2y D«y 



dx 



o L r(a + i) 

when subject to the boundary conditions 

2/(0) = 0, y(l) = l. 

Here, y := x a ,x € [0, 1]. The fractional Euler-Lagrange associated to this problem is 



x 1 - 



T(a + V, 



( D° y y-2y D«y 



2x a 



r(a + i) 



oD%y - 2y 



0. 



(20) 



Replacing y by y, and since §D%y = T(a + 1), we conclude that y is a solution of (f20l) . 

Consider now the following problem: what is the order of the derivative a, such that ^(y) 
attains a maximum or a minimum? In other words, find the extremizers for ip(a) = ^f(y). 
Direct computations show that 



Jo 



^(a) = / [x a T(a + l)Y dx. 

Evaluating its derivative, 
rl a 

if/ ( a ) = / -r- [x a T{a + l)] 2 dx 
Jo da 



2x a T(a + 1) 



/>oo 

x a In x r(a + l)+x a t a In t e^dt 
Jo 



dx. 



We have that a 0.901 is a solution of the equation ip'(a) = 0, and such value is precisely 
where ^?(y) attains a minimum. 
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0.25- 



0.0- 



I I I I | I I I I | I I I I | I I I I | 
0.0 0.25 0.5 0.75 1.0 

alpha 



Figure 2: Graph of \&(y) for a G [0, 1] 



More generally, consider the functional 

*(»,«)= / L(x,y(x), a D%y(x))dx. (21) 

J a 

Functional ()21 j) contains the left Riemann-Liouville derivative only, but we can consider func- 
tional containing right Riemann-Liouville derivatives or both in a similar way. Let h be a 
curve such that h(a) = hip) = 0, 5 be a real number, and (y,a) be an extremal for <£. Then, 

$(y + /i,a + <T) - $(y,a) 

fb aj ar ar 

= J ^■h+ — . a D° +5 h + — -( a D«+ 5 y- a D«y)dx + 0\(h,6)\ 2 . 

For 5 = 0, using the fractional integration by parts formula and the fundamental lemma 
of the calculus of variations, we obtain the known fractional Euler-Lagrange equation: 

g-(x,y(x), a D«y(x)) + x DZ — (x,y(x), a D%y(x)) = 0. 



For h = 0, we obtain the relation 

du 



x,y(x), a D%y{x))<f>' {a)dx = 0, 



where 4>{a) = a D x y(x). In summary, we have: 

Theorem 12. If(y,a) is an extremal of § given by (|2ip . satisfying the boundary conditions 
y(a) = and y(b) = y&, then y satisfies the system 



( BT 81 

— (x,y(x), a D«y(x))+ x DZ — (x,y(x), a D a x y(x)) = 

— (x,y(x), a D%y(x))(jj \a)dx = 

^ Ja OU 



(22) 



where <f>{a) = a D x y(x). 
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In the previous example the solution obtained satisfies system (|22p since 



— (x,y(x), a D%y(x))) = Q. 
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